高雷诺数下非混相Rayleigh-Taylor不稳定性的格子Boltzmann方法模拟

本文采用相场格子Boltzmann方法研究了竖直微通道内中等Atwoods数流体的单模Rayleigh-Taylor不稳定性问题,系统分析了雷诺数对相界面动力学行为以及扰动在各发展阶段演化规律的影响.数值结果表明高雷诺数条件下,不稳定性界面扰动的增长经历了四个不同的发展阶段,包括线性增长阶段、饱和速度阶段、重加速阶段及混沌混合阶段.在线性增长阶段,我们计算获得的气泡与尖钉振幅符合线性稳定性理论,并且线性增长率随着雷诺数的增加而增大.在第二个阶段,我们观察到气泡与尖钉将以恒定的速度增长,获得的尖钉饱和速度略高于Goncharov经典势能模型的解析解[Phys.Rev.Lett.200288134502],这归因于系统中产生了多个尺度的旋涡,而涡之间的相互作用促进了尖钉的增长.随着横向速度和纵向速度的差异扩大,气泡和尖钉界面演化诱导产生的Kelvin–Helmholtz不稳定性逐渐增强,从而流体混合区域出现许多不同层次的涡结构,加速了气泡与尖钉振幅的演化速度,并在演化后期阶段,导致界面发生多层次卷起、剧烈变形、混沌破裂等行为,最终形成了非常复杂的拓扑结构.此外,我们还统计了演化后期气泡与尖钉的无量纲加速度,发现气泡和尖钉的振幅在后期呈现二次增长规律,其增长率系数分别为0.045与0.233.而在低雷诺条件下,重流体在不稳定性后期以尖钉的形式向下运动而轻流体以气泡的形式向上升起.在整个演化过程中,界面变得足够光滑,气泡与尖钉在后期的演化速度接近于常数,未观察到后期的重加速与混沌混合阶段.     

对流扩散方程的缩减基有限元法的后验误差

将缩减基(RB)方法和有限元方法相结合,在保证偏微分方程的有限元离散格式具有足够高精确度前提下,能够大幅度地降低有限元离散格式的维数,从而大大降低计算中内存容量和计算时间的消耗.针对对流扩散方程建立基于RB方法的Crank-Nicolson有限元离散格式,并给出后验误差估计结果.     

线性方程组求解与向量组规范化的初等变换方法

本通过对线性方程组的系数矩阵的行与列的初等变换给出了求解线性方程组的方法，并通过对矩阵的初等变换给出了向量组正化的方法。     

数字散斑相关方法及其在碳纤维复合材料压力容器变形测量中的应用

提出了一种新的基于图像灰度梯度迭代的数字散斑相关方法(DSCM,digital speckle correlation method)。通过使用DSCM测量碳纤维复合材料压力容器在水压下的局部区域的位移场和应变场,分析了复合材料压力容器的轴向和环向的变形特征,为碳纤维复合材料压力容器的优化设计提供了理论和实验依据。     

对问题1627的质疑

1.对于问题1626的繁琐解答,我们陆续收到了很多读者的来信,对其予以改进.由于来函较多,而此问题的改进解法确又非常简单(通过连线,利用勾股定理或余弦定理即可解出.当然还有其它简单解法),在此我们就不再选用新的解法刊出了.另外,北京的贺信淳老师对我们所选题目也进行了批评,我们诚恳地表示接受并对贺老师及所有给我们来信的老师表示衷心的感谢.今后,我们会继续努力,加强自身建设,争取和作者们一道把问题栏越办越好.真诚地欢迎关心爱护我们的读、作者随时给我们提出宝贵的意见及建议.感谢大家对问题栏的支持.2.对于问题1627解答的错误,读者们给出了很多不同的正确解答.根据来稿时间及解答方法,我们选用了湖南师大梁红梅等三位老师合写的文章,对原解答予以纠正.对其它提供正确解答的老师,我们在此一并表示感谢.     

一类非线性Schr(o)dinger方程的守恒差分法与Fourier谱方法

考察了一类带导数项的非线性Schrodinger方程的周期边值问题，提出了一种守恒的差分格式，在空间方向上采用Fourier谱方法，证明了格式的稳定性和收敛性．数值试验得到了与理论分析一致的结果．     

氦原子与非对称替代氢分子碰撞(E=0.3eV)角分布计算

用密耦方法计算了非对称同位素替代分子与氦原子碰撞(He-HD,HT,DT)转动激发,当入射能量E=0.3 eV时,得到了上述碰撞体系的微分散射截面或角分布.基于上述理论计算,讨论了原子与双原子分子碰撞的同位素效应.     

理论研究C8H80-（H2O）n（n=1～5）团簇

应用密度泛函B3LYP／6—31＋G（d，p）方法对C8H80-（H2O）n（n=1～5）团簇这种弱相互作用体系进行垒自由度能量梯度优化，得到该系列团簇的稳定蛄构．计算结果表明。在该系列二元团簇中，一方面水分子数目的多少对苯基丙酮分子的结构影响很小，另一方面由于苯基丙酮分子的存在，破坏了团簇中水分子的对称性结构，在团簇内部极力形成O—H—O这样弯曲的有方向性的氢键．对苯基丙酮-水这样结构复杂的团簇，指认光谱的难度非常大，本文只讨论了与C=O有关的振动峰和水分子的对称伸缩振动的最强峰．